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1 Problem 

Deduce the form of a cylindrically Symmetrie plane electromagnetic wave that propagates in 
vaeuum. 

A scalar, azimuthally Symmetrie wave of frequency uj that propagates in the positive z 
direction could be written as 

^(r,t) = /(p)e^-^, (1) 



where p = \Jx 2 + y 2 . Then, the problem is to deduce the form of the radial funetion f(p) 
and any relevant condition on the wave number k z , and to relate that scalar wave funetion 
to a complete Solution of Maxwell's equations. 

The waveform ([]]) has both wave velocity and group velocity equal to uo/k z . Comment 
on the apparent superluminal character of the wave in case that k z < k = u/c, where c is 
the speed of light. 

2 Solution 

As the desired Solution for the radial wave funetion proves to be a Bessel funetion, the 
cylindrical plane waves have come to be called Bessel beams, following their introduetion 
by Durnin et al. [[1], 0]. The question of superluminal behavior of Bessel beams has recently 
been raised by Mugnai et al. ||. 

Bessel beams are a realization of super-gain antennas j|, ^|, |5| in the optical domain. A 
simple experiment to generate Bessel beams is described in j7|. 

Sections 2.1 and 2.2 present two methods of Solution for Bessel beams that satisfy the 
Helmholtz wave equation. The issue of group and signal velocity for these waves is discussed 
in sec. 2.3. Forms of Bessel beams that satisfy Maxwell's equations are given in sec. 2.4. 

2.1 Solution via the Wave Equation 

On substituting the form ([!]) into the wave equation, 

< 2 > 

we obtain 

^ + -f + (k 2 -kl)f = 0. (3) 
ap z p dp 

This is the differential equation for Bessel funetions of order 0, so that 



f{p) = MKp), (4) 
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where 

K + kl = k 2 . (5) 
The form of eq. (||) suggests that we introduce a (real) parameter a such that 

k p = ksina, and k z = kcosa. (6) 
Then, the desired cylindrical plane wave has the form 

iß(r,t) = J {ksinap)e i{ - kcoaaz - ut \ (7) 

which is commonly called a Bessel beam. The physical significance of parameter a, and that 
of the group velocity 

du üü c 
9 dk z k z p cosa 

will be discussed in sec. 2.3. 

While eq. (|?p is a Solution of the Helmholtz wave equation (|2|), assigning ^(r, t) to be a 
single component of an electric field, say E x , does not provide a füll Solution to Maxwell's 
equations. For example, if E = then V • E = dip/dx ^ 0. Bessel beams that satisfy 
Maxwell's equations are given in sec. 2.4. 

2.2 Solution via Scalar Diffraction Theory 

The Bessel beam ([?]) has large amplitude only for \p\ < 1/fcsina, and maintains the same 
radial profile over arbitrarily large propagation distance z. This behavior appears to con- 
tradict the usual lore that a beam of minimum transverse extent a diffracts to fill a cone of 
angle 1/a. Therefore, the Bessel beam (0) has been called "diffraction free" ||. 

Here, we show that the Bessel beam does obey the formal laws of diffraction, and can be 
deduced from scalar diffraction theory. 

According to that theory ||, a cylindrically Symmetrie wave f(p) of frequency u at the 
plane z = propagates to point r with amplitude 

U ff p i(kR—u>t) 

1>(r, t ) = ^-jj P'dp'dHip')-^—, (9) 

where R is the distance between the source and Observation point. Defming the Observation 
point to be (p, 0, z), we have 

R 2 = z 2 + p 2 + p 2 - 2pp'cos0, (10) 

so that for large z, 

2z v ; 

In the present case, we desire the amplitude to have form ([[]). As usual, we approximate 
R by z in the denominator of eq. (|9]), while using approximation ( pÄ|) in the exponential 
factor. This leads to the integral equation 

U ikz ikp 2 /2z roo , r 2n 

f(p)e ik * z = —-— / p'dp'f(p')e ikp /2z / dct)e- ikpp co ^ /z 



2ixi z 

gikZgikp 2 /2z „oo 







p'dp'f(p')J (kpp'/z)e^ 2 / 2z , (12) 
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using a well-known integral representation of the Bessel function J . 

It is now plausible that the desired eigenfunction f(p) is a Bessel function, say J (k p p), 
and on Consulting a table of integrals of Bessel functions we find an appropriate relation |J, 

p'dp'J (k pP ')J (kpp'/z)e ikp ' 2/2z = ^e- ikp2/2z e- ik > /2k J (k p p). (13) 

Comparing this with eq. (0), we see that f(p) = Jo(k p p) is indeed an eigenfunction provided 
that 2 

k z = k-^. (14) 



2k 



Thus, if we write k p = fcsina, then for small a, 



k z « fc(l - a 2 /2) ps fccosa, (15) 

and the desired cylindrical wave again has form (p]). 

Strictly speaking, the scalar diffraction theory reproduces the "exact" result (|7]) only for 
small ol. But the scalar diffraction theory is only an approximation, and we predict with 
confidence that an "exact" diffraction theory would lead to the form fl7|) for all values of 
parameter ol. That is, "diffraction- free" beams are predicted within diffraction theory. 

It remains that the theory of diffraction predicts that an infinite aperture is needed to 
produce a beam whose transverse profile is invariant with longitudinal distance. That a 
Bessel beam is no exception to this rule is reviewed in sec. 2.3. 

The results of this section were inspired by [|IIJ. One of the first Solutions for Gaussian 



laser beams was based on scalar diffraction theory cast as an eigenfunction problem [11 



2.3 Superluminal Behavior 

In general, the group velocity (§) of a Bessel beam exceeds the speed of light. However, this 
apparently superluminal behavior cannot be used to transmit signals faster than lightspeed. 

An important step towards understanding this comes from the Interpretation of parameter 
a as the angle with respect to the z axis of the wave vectors of an infinite set of ordinary 



plane waves whose superposition yields the Bessel beam | 12| . To see this, we invoke the 



integral representation of the Bessel function to write eq. ([7]) as 
ip(r,t) = J {ksmap)e i{kcoBaz - ut) 



1 r 2ir 

"06 ■ s ' n a x cos s ^ na v s ™ < j^+' c cos a z—u>tj (16) 



2tt Jo 

1 r 2n 
= — / rf0e i(q - r ^^ ) , 
2ix Jo 

where the wave vector q, given by 

q = k(sin ol cos0, sin a sin 0, cosa), (17) 
makes angle a to the z axis as claimed. 



3 



We now see that a Bessel beam is rather simple to produce in principle 0. Just superpose 
all possible plane waves with equal amplitude and a common phase that make angle a to 
the z axis, 

According to this prescription, we expect the z axis to be uniformly illuminated by the 
Bessel beam. If that beam is created at the plane z — 0, then any annulus of equal radial 
extent in that plane must project equal power into the beam. For large p this is readily 
confirmed by noting that Jq (k sin ap) ~ cos 2 (k sin a p + 5)/(k sin a p) , so the integral of the 
power over an annulus of one radial period, Ap = 7ij(ksina), is independent of radius. 

Thus, from an energy perspective a Bessel beam is not confined to a finite region about 
the z axis. If the beam is to propagate a distance z from the plane z = 0, it must have radial 
extent of at least p = ztana at z — 0. An arbitrarily large initial aperture, and arbitrarily 
large power, is required to generate a Bessel beam that retains its "diffraction-free" character 
over an arbitrarily large distance. 

Each of the plane waves that makes up the Bessel beam propagates with velocity c along 
a ray that makes angle a to the z axis. The intersection of the z axis and a plane of constant 
phase of any of these wave moves forward with superluminal speed cj cosa, which is equal 
to the phase and group velocities 

This superluminal behavior does not represent any violation of special relativity, but is 
an example of the "scissors paradox" that the point of contact of a pair of scissors could 
move faster than the speed of light while the tips of the blades are moving together at 
sublightspeed. A ray of sunlight that makes angle a to the surface of the Earth similarly 
leads to a superluminal velocity cj cosa of the point of contact of a wave front with the 
Earth. 

However, we immediately see that a Bessel beam could not be used to send a signal from, 
say, the origin, (0,0,0), to a point (0,0, z) at a speed faster than light. A Bessel beam at 
(0, 0, z) is made of rays of plane waves that intersect the plane z = at radius p = ztana. 
Hence, to deliver a message from (0, 0, 0) to (0, 0, z) via a Bessel beam, the Information must 
first propagate from the origin out to at least radius p = ztana at z = to set up the beam. 
Then, the rays must propagate distance zj cosa to reach point z with the message. The 
total distance traveled by the Information is thus z(l + sina)/ cosa, and the signal velocity 
v s is given by 

(18) 



1 + sin a 

which is always less than c. The group velocity and signal velocity for a Bessel beam are 
very different. Rather than being a superluminal carrier of Information at its group velocity 
cj cosa, a modulated Bessel beam could be used to deliver messages only at speeds well 
below that of light. 

2.4 Solution via the Vector Potential 

To deduce all components of the electric and magnetic fields of a Bessel beam that satisfies 
Maxwell's equation starting from a Single scalar wave function, we follow the Suggestion of 
Davis [13| and seek Solutions for a vector potential A that has only a Single component. We 
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work in the Lorentz gauge (and Gaussian units), so that the scalar potential $ is related by 

(19) 



. 1<9$ ft 
V-A+-— = 0. 

c öt 



The vector potential can therefore have a nonzero divergence, which permits Solutions having 
only a single component. Of course, the electric and magnetic fields can be deduced from 
the potentials via 

1 BA 

E = -V$--^, (20) 



and 



c dt 



B = V x A. 



(21) 

For this, the scalar potential must first be deduced from the vector potential using the 
Lorentz condition (|T9|). We consider waves of frequency uj and time dependence of the form 



-iujt 



so that d<&/dt = —ikQ. Then, the Lorentz condition yields 



k 



-V- A, 



and the electric field is given by 



E = ik 



A + -V(V-A) 



(22) 



(23) 



Then, V • E = since V 2 (V • A) + k 2 (V ■ A) = for a vector potential A of frequency to 
that satifies the wave equation (g), etc. 

We already have a scalar Solution (^) to the wave equation, which we now interpret as 
the only nonzero component, Aj, of the vector potential for a Bessel beam that propagates 
in the +z direction, 



Aj(r,t) = ip(r,t) oc J (k sin a p)e 



i(k COS OL Z — Lüt) 



(24) 



We consider five choices for the meaning of index j, namely x, y, z, p, and 0, which 
lead to five types of Bessel beams. Of these, only the case of j — z corresponds to physical, 
azimuthally Symmetrie fields, and so perhaps should be called the Bessel beam. 



2.4.1 j = x 

In this case, 



V- A 



di/j 
dx 



k sin ol x 
P 



J\{k sin a p)e 



i(k COS OL Z—Lüt) 



(25) 



In calculating V(V • A) we use the identity J[ = (J — Jz)/2. Also, we divide E and B by 
the factor ik to present the results in a simpler form. We find, 



Er 



E- 



y (g)- 

sin 2 ol xy 



sin 2 a 



fr 



Q 2 



i(k cos a z—ut) 



J -^-\{Mq)-Uq)) 



3 i(fe COS OL Z — Lüt) 



(26) 



x 



-i sm2a—J 1 (g)e i(kcosaz - Lüt \ 
2p 



where 

q = k sin a p, (27) 

and 

B x = 0, 

B y = cos et Jo^e^ 08 02 "^, (28) 
B z = -i sin a^J^gy^™ 3 «^. 

A Bessel beam that obeys Maxwell's equations and has purely x polarization of its electric 
field on the z axis includes nonzero y and z polarization at points off that axis, and does not 
exhibit the azimuthal symmetry of the underlying vector potential. 

2.4.2 j = y 

This case is very similar to that of j = x. 

2.4.3 j = z 

In this case the electric and magnet fields retain azimuthal symmetry, so that it is convenient 
to display the p, <fi and z components of the fields. First, 

V- A = ^ = ^cosaJo(fcsinap)e i(fccosaz -^ ) . (29) 

Then, we divide the electric and magnetic fields by k sin a to find the relatively simple forms: 

E p = cos« J 1 (ß)e <(fcC0Ba *- w * ) , 

= 0, (30) 
E z = isiYiaJ {Q)e i{kcosaz -^\ 

and 

B P = 0, 

B^ = Ji(e)e < ( fc0 ° saÄ - wt \ (31) 
B z = 0. 

This Bessel beam is a transverse magnetic (TM) wave. The radial electric field E p vanishes 
on the z axis (as it must if that axis is charge free), while the longitudinal electric field E z 
is maximal there. Cylindrically Symmetrie waves with radial electric polarization are often 
called axicon beams |14] . 



2.4.4 j=p 

In this case, 

1 dpip 



V- A 



p dp 



J (k sin ap) 

hina J\ [k sm a p) 



g*(fc COS OL Z — U)t) (32) 
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After dividing by ik, the electric and magnetic fields are 



and 



E p = \ J (q) - sin 2 a 



j,(k COS OL z—wt) 



E<f> 
E- 



i cos a sm a 



Mo) 



Mo) 



3 i(fe COS OL Z—Lüt) 



(33) 



B p = 0, 

= cosaJ {g)e i(kcoaaz - ut \ 
B z = 0. 

The radial electric field diverges as 1/p 2 for small p, so this case is unphysical. 



(34) 



2.4.5 j = <f> 

Here, 

After dividing by ik, the electric and magnetic fields are 



E p = 0, 

E<ß = J {ß)e^ kcosaz - Wt \ 
E t = 0, 



(35) 



(36) 



and 



B^ = 
B ? = 



cosaJ {o)e l{kcosaz ~ u)t) 

Mo) 



% sm a. 



Q 



Ji(ß) 



i(k COS OL Z — Lüt) 



(37) 



These fields are unphysical due to the finite value of at p = 0, and the divergence of B z 
as p — > 0. 
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